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We revisit the canonical Rayleigh-Taylor instability and investigate the case of a 
thin film of fluid upon the underside of an inclined plane. The presence of a natural 
flow along the plane competes with the conventional droplet forming instability. In 
particular, experiments reveal that no drops form for inclinations greater than a crit¬ 
ical value. These features are rationalized in the context of the absolute/convective 
analysis conducted in this article. 


I. INTRODUCTION 

The Rayleigh-Taylor instability (RTI) occurs when a fluid is momentarily sustained above 
a fluid of lesser density 1 ! Examples are found in a variety of situations ranging from the 
simple case of water sitting on oil in a kitchen glassware to supernova explosions 2 . This 
interfacial instability has been thoroughly studied and has served as a foundation for tackling 
fundamental issues such as the break up of free-surface flowP (drop formation) and pattern 
generation 4 . Applicative components of this instability also arise due to its natural tendency 
to enhance the mixing of liquid specie^! Additionally, the RTI is a prime concern when 
coating surfaces, be it with paint or a lubricant, as the instability may lead to coating 
irregularities or to the detachment of droplets for thick coatings. As such, many studies 






FIG. 1. (a) A thin film of initial thickness hi flows on the underside of a transparent inclined plane. 

The angle a denotes the inclination of the substrate with respect to the gravity g. (b) A typical 
experimental observation for hi — 0.52 I c and a = 81°. (c) No drops are formed for hi — 0.55 I c 
and a = 60.5°. (d) For large enough values of a , droplets form according to the Rayleigh-Taylor 
instability. The black arrows indicate the overall direction of the flow. 
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have concentrated on means of controlling or suppressing the growth of pendant drops when 
applying temperature gradients, electric fields, or surface tension gradients (see Ref . 5 for 
example). In the following, we revisit the iconic problem of a thin liquid film spread upon 
the underside of a flat surface and investigate the effect of the substrate orientation on the 
instability. 

In the canonical RTI of a thin viscous fluid layer coated on the underside of a horizontal 
surface, the interface perturbations are found to grow exponentially in time and generate 
drops arranged in various patterns 4 . The instability development is set by the competition 
of surface tension effects and gravitatio nal eff ects such that the distance between droplets 
is found to be A = 2i\\[2£ c where £ c = ^/j/pg is the capillary length of the problem and 7 , 
p, and g denote the thin layer surface tension, its density, and the acceleration of gravity, 
respectively. This natural length scale is well recovered with linear arguments since A is the 
wavelength of the most unstable mode predicted by linear stability analysis. 

Herein, we extend the study to cases where the surface is tilted and set at an angle a ^ 7 r /2 
from the vertical (see Fig. [l] a = 7 r /2 corresponds to the horizontal case). In particular, 
we focus on how the main features of the instability are affected by the base flow along the 
substrate. For values of a ~ n/2 we observe that the interface deforms and eventually forms 
droplets according to the RTI (Fig. [ljb,d)). For values of a < 7 r/ 2 , deformations of the 
surface are greatest in the downstream direction where the fluid accumulates; this is also 
the area where the first droplets form. Soon after, more droplets form and drip from the 
upward area as well. However, building upon the known characteristics of flows on vertical 
walls - be it from the application of varnish on furniture in daily experience or else when 
referring to the literature on plates extracted from liquid baths 6 - one may anticipate that 
the inclination of the substrate is likely to affect the instability to a point where it could 
disappear, say at a critical angle a c . This angle a c has been observed experimentally (see 
Fig. 0 c)) as explained in the results Sec. |IV| We will see that traditional temporal linear 
stability arguments, which are sufficient in the horizontal case, fail to predict the value 
of a c . The experiments are instead rationalized by addressing the absolute or convective 
nature of the aforementioned linear stability results as detailed hereafter. Note that the 
Kapitza instability^ is here neglected owing to the low Reynolds numbers characterizing 
the experiment. Hence, we solely focus our analysis on the RTI. 


II. THE FLOW UNDER AN INCLINED PLANE 



FIG. 2. A thin film of initial thickness hi flows under an inclined plane in the direction e x . The 
angle a denotes the inclination of the substrate with respect to the direction of gravity g. 
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We consider the flow of a thin fluid film of initial thickness hi placed underneath an 
inclined plane (Fig. [l]). The film, of density p, viscosity p, and surface tension 7 is subject 
to the acceleration of gravity g. Of particular importance herein is the angle a > 0 between 
the substrate and the vertical. The thickness hi is supposed smaller than the typical length 
scale over which the coating spreads, i.e., the problem is assumed invariant in the flow 
direction. Additionally, mass conservation indicates that the flow velocity u = (u,v) is 
essentially one dimensional (v ~ u x h/L <C u). Assuming low Reynolds number, one may 
write the lubrication equations 9 for this flow yielding: 

. N pg ( . dh 7 8k \ , 

u(x,y) = — l-sma- - — — - cosa 1 y{y - 2 h) (1) 


where k ~ is the interface curvature. Integrating eq. |l]) leads to finding the flow rate 
Q oc h 3 , which combined to mass conservation yields: 
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The trivial solution h = hi satisfies eq. ([ 2 ]) and is valid for any value of hi <C L. In that 
case, the Poiseuille like flow in the x-direction is such that the free surface velocity writes: 


Ui = cos a — h 2 

p 


( 3 ) 


We now turn to examining the stability of the uniform solution to linear perturbations. 
Assuming a normal decomposition h(x, t ) = hi+ehd(x, t) with e< 1 and hd(x, t ) oc 
(see Fig. [ 2 ]), one obtains the following dispersion relation for uj from eq. 
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The real part of uo represents the advection of the wave and is found to be equal to rq, 
the velocity at the interface of a uniform film. The imaginary part of cj, denoted cr, is the 
temporal growth rate of the instability, whose sign sets the stability of the perturbation. 
The flow is either stable, marginally stable, or unstable for a < 0, a = 0, and a > 0 
respectively. 

As one may observe in Fig. [3j the flow is always found to be linearly temporally unstable 
for strictly positive values of a and the most unstable modes are: 


k 


max 


1 


Vsum. 


( 5 ) 


They correspond to a wave length A max = 27rv / 2^ c /v / sum and comments may be made 
regarding this expression. First, when substituting a = n/2 in eq. one recovers the 
wavelength of the RTI of a thin film underneath an horizontal plate, X a=n /2 = 27r\/2 l c , 
which sets the distance between neighboring droplets (on the order of 2cm here). Second, 
with respect to a is limited when a ~ tt/2 (see the y/sin a in eq 


the sensitivity of A 

and the proximity of the higher most two curves in Fig.J^a)). Finally, a crude interpretation 
of eqs. 0 and © would lead to the conclusion that only perfectly vertical walls prevent 
drops from forming. In Fig. |3ja) the dot-dash curve corresponding to a = 0 is the only 
one fully enclosed in the stable region. In other words, the value of the earlier defined 
critical angle would be a c = 0. This is in contradiction with our experimental observation 
presented in Fig. [ljc) where a ~ 7r/3 a nd n o drops form. This paradox will be solved by 


the absolute/convective analysis in Sec. Ill but a first insight in the mechanism at stake 


may be gained when evaluating the typical time scales of the problem. 
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(a) 


(b) 



FIG. 3. (a) Growth rates for inclinations varying from 7r/2 to 0 with a 7r/20 increment. The most 
unstable mode is marked by a black circle. All plain lines denotes unstable cases and the dotted 
dash line (corresponding to a = 0) is found marginally stable. For a = 7r/2 the maximum growth 
rate corresponds to k max = ic 1 /V 2. (b) Characteristic times of the instability growth (r max , plain 
line) and of its advection (r/). The dash line corresponds to hi/£ c = 0.5. 


We now proceed to evaluate and compare the characteristic times of the instability and 
of the flow. Let us define r max to be the inverse of the growth rate of the most unstable 
mode, i.e., r max = a(k ma , x )~ 1 . Using eqs. 0 and ([H]) we find that: 

Tnax CK ~ o (6) 

sin a 

Similarly, Tf is defined as the time of advection over the wavelength A max , i.e., Tf ~ A max /iq 
yielding: 

V K ( 7 ) 

cosaysma 

As shown in Fig. [3|b), Tf and r max are respectively increasing and decreasing with a. 
In particular, they intersect for a given angle denoted a*. The practical implications are 
as such: for small angles a < a* we anticipate that the instability will be dominated by 
the flow since the time required for the instability to form (r max ) is much larger than the 
typical time needed for the same instability to flow down Tf. Conversely, at large angles the 
instability develops much faster than the fluid flows, i.e., we anticipate that drops will form 
as if there were no flow (in the limit of very large angles a ~ 7r/2). The distinction between 
these two regimes may be made rigorously in the context of a spatio-temporal analysis as 
detailed in the following. 


III. ABSOLUTE AND CONVECTIVE INSTABILITIES 

In this section, we make use of the absolute/convective instability concepts to analyze 
the RTI on the underside of an inclined surface. The distinction between absolute and 
convective instability has been pioneered in plasma instabilities 10 11 and then more recently 
applied to fluid-dynamical instabilities in parallel and weakly spatially developing open 
fk>wd l2 l The nature of a given flow depends on the large-time asymptotic behavior of the 
linear impulse response; the flow is convectively unstable if the amplified disturbances move 
away from the source, conversely, the flow is absolutely unstable when amplified perturba¬ 
tions invade the entire flow. In contrast to the temporal stability problem where the axial 
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wavenumber k is real and one seeks a complex frequency uj, the absolute/convective nature 
of the instability is determined by applying the Bri ggs-Be rs zero-group velocity criterion to 
the dispersion relation for fully complex (fc, uj) pairs 10 11 . In order to determine the transi¬ 
tion from convective to absolute instability, it is sufficient to detect the occurrence of saddle 
points in the characteristics of spatio-temporal instability waves, i.e., a complex value of 
the wavenumber fco such that = 0. The discrimination between absolutely and 

convectively unstable flows was shown to play a crucial role in accounting for the occurrence 
of synchronized self-sustained oscillations in a variety of spatially developing shear flows, 
such as single phase wakes, hot jets, and counter-flow mixing layers 1 ^. But it also explains, 
for instance, the transition from dripping to jetting in two-phase immiscible microfluidic 
co-axial injectors 1 ^, as well as different regimes occurring in the pearl forming instability of 
a film flowing down a fibeiW 

In anticipation of the following analysis we rewrite eq. © in a convenient dimensionless 
form: 

uj = i [k 2 — fc 4 ^ + uk , (8) 

where the wavenumber and angular frequency have been non-dimensionalised with the char¬ 
acteristic space and time scales as: 


k = k- 
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A/sum 

The non-dimensional interface velocity reads: 


and uj = uj 
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( 9 ) 


u = 


Vsum tan a hi 


( 10 ) 


Hereafter, both the wavenumber k and the wave-function uj are complex numbers. De¬ 
noting k = k r + iki and injecting this form in eq. ^ yields: 

LO r = k r (-4fcf + ki - 2) + Uj (11) 

U>i = —kf + k^ ^6kr — 1 ) + kill — kr + kr (12) 

The complex group velocity writes v g = and we denote ujq and ^0 the absolute fre¬ 
quency and wavenumber, respectively. They are defined bv the zero-group-velocity condi¬ 
tion Vg(ko ) = 0, such that uj$ = uj(ko). Differentiating (I12| with respect of k r and ki and 
canceling both equations, one obtains: 


u = -32/co,i - 4/c 0 , i 

^0,1 = 1/4 - 2^o 5 i - 24 ^o 5 • 


(13) 

(14) 


Note that the values of fco,* an( 4 ,* ai *e directly set by the value of the parameter u , itself 
function of the physical parameters hi/t c and a (see eq. (10)). In turn, the value of 
relative to zero sets the nature of the flow. When ujq^ < 0 the flow is convectively unstable 
and conversely, it is absolutely unstable when a) 0 ,* > 0- I n both cases the flow is linearly 
unstable, as found earlier in Sec. |II[ However, those two cases lead to completely different 
scenarios, which are the main focus of the rest of the article. We now search for the roots of 
eq. (14), which combined with eq. ©, yields to the exact expression of the critical value 
of u = u* leading to a transition between the convective and absolute regimes, namely: 


u = 


V27 


34 + UV7 ~ 1.62208. 


(15) 
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FIG. 4. Spatio-temporal evolution of two wave packets generated by a Dirac perturbation in x — 0 
at t — 0 + for hi/i c ~ 0.67 and (a) a — 7r/4 and (b) a = 7r/2.2. Both cases are unstable, but 
(a) is convectively unstable, the perturbation decreases with time in x — 0, and conversely, (b) 
is absolutely unstable as the perturbation grows exponentially in the laboratory frame. Dashed 
lines correspond to the fronts of the perturbed wedges, (c) shows the time evolution of the relative 
amplitude of the perturbation at the origin in both cases. 


With a dispersion relation taking the form found in eq. (J8J) , the leading v j and receding 
vj front velocities of the wedge of the impulse response can be directly determined as the 
sum of the inherent advection u and db u* : 


_± ~~, Av f 

VJ = u±u = u± —- jL 
J 2 


(16) 


where A Vf = v^ — vj = 2 u* is the width of the perturbed region. This concept is best 
illustrated when following the time evolution of a wave packet generated with a localized 
initial perturbation. In Fig. [4j we report the time evolution of two packets generated by a 
Dirac impulse in x = 0 at t = 0 and observed in the laboratory frame. First, a global view 
of the flow as a function of time is presented in the two panels in Fig. |4ja) and (b). They 
correspond to Co o,i < 0 and uj o,i > 0, respectively. Second, the time evolution at the impulse 
location, ft,(0, £), is presented in Fig. [4^c). When < 0, the perturbation increases with 
time but is carried away faster than it spreads across the flow and the upstream front velocity 
vj of the perturbed wedge is positive. As a consequence, the perturbation at the source 
point (x = 0) decreases with time such that h(0,t) is quickly indistinguishable from the 
vertical axis, and the instability is convective. Conversely, if d>o,i > 0, the perturbation grows 
exponentially in the laboratory frame despite the presence of the flow. The upstream front 
velocity vj of the perturbed wedge is negative, so the wave-packet counter-propagates to 
invade the all space, and the instability is absolute (Figure [ 4 ) 3 ) . These examples emphasize 
that the absolute/convective transition is found when vj = 0, i.e., u = u*. 

Recalling from eqs. 


stant, one may derive t 
transition: 


lOj) and ( 15]) that ft is a decreasing function of a and that u* 
re critica 


value of the inclination a* yielding the absolute/convective 


Vsum* tancC = 


34 

u*hi 


(17) 


The existence of such an angle had been anticipated in the time scale analysis in Sec. |Tl 
using temporal arguments. However, the presented spatio-temporal analysis constitutes a 
rigorous derivation of its exact expression. Note that these results are consistent with the 
mechanism depicted in Sec. [TT} Representing eq. (16) with dimensional terms, the existence 
of a critical angle results from the competition of the advection velocity iq = cos 
which increases with an increasing inclination (decreasing values of <a), and the wedge 
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spreading of the unstable region: 


A v f 


= 2 u 


~*Pghl sin 3 / 2 a 


34 M 


(18) 


which decreases with an increasing inclination (decreasing values of a). A critical angle is 
unavoidable as an increasing inclination of the plane therefore simultaneously increases the 
advection and decreases the growth rate. 


We turn back to examining the critical angle in eq. ( |17| ). This implicit expression is 
solely a function of the ratio hi/t c , such that surface tension and gravity effects influence 
the problem through the capillary length £ c . On the other hand, the fluid viscosity /i does 
not enter in eq. (17). This could have been anticipated as fi intervenes linearly both in the 
typical flow rate and in the typical droplet formation time, such that it does not affect their 
relative values. A practical consequence of eq. © is the state diagram shown in Fig. |] 



FIG. 5. State diagram of the flow under an inclined plane in the parameter space (hi/£ Cl a). Plotted 
in black are the values of the critical angle a* delimiting the absolute and convective domains. The 
two stars indicate the physical parameters used in Fig. [4] 


where the plot of a* as a function of hi/£ c may be found. For a given value of hi /£ c , the 
flow is absolutely unstable for angles a > a* and convectively unstable for a < a*. Note 
that for a = 7t/ 2, all thin films with non-zero thickness are absolutely unstable since such 
an angle forbids any flow parallel to the substrate to compete with the instability. 

In summary, the flow under an inclined plane is always linearly temporally unstable 
(unless a = 0). However, two cases are possible. Depending on the respective values 
of a and hi/£ c the flow is either convectively or absolutely unstable. In the first case, 
perturbations are found to decrease in the laboratory frame, despite the unstable nature 
of the flow, as they are carried away by advection. In the second, the instability is strong 
enough to overcome the flow. We now examine experimentally the consequences of such 


scenarios. 






FIG. 6. (a) Experimental results. Each experiment is represented by the triplet ( hi/i c , a , N) and is 
color coded according to the scaled number of droplets formed N. (b) Experimental data compared 
to the theoretical expression a *. 


IV. EXPERIMENTS 

The experiments were conducted on a test stand which consists of a base plate, two 
columns, and two mounting plates where the clamp assemblies that hold the experimental 
surface in place are free to rotate (Fig. [lja)). The working fluid used in the experiments is 
Castor Oil (Hanseller AG). The value of the capillary length £ c = 1.91mm was measured 
using Axisymmetric Drop Shape Analyst 5 on a pendant drop experiment. Plexiglas was 
chosen to serve as support surfaces for the experiments. A grid was included in the setup 
to enhance the distortions of the interface. The grid was aligned parallel with the surface 
and placed far below the flow so as not to induce any perturbations. 

A typical experiment is divided in two steps: preparation and testing. First, a fluid film is 
prepared on the solid surface: the surface is made horizontal and a mass m of oil is deposited 
in the center of the surface so that it will spread over time. Its area A increases with time, 
A oc £ 1//8 , while its thickness remains uniform apart from the boundaries 16 . Second, at a 
given time t — ti the surface is gently rotated into the testing position where the fluid is 
beneath the surface, which makes an angle a with the vertical axis e z (Fig. [lja)). At the 
beginning of the testing, the thickness of the film writes: hi = ^ A ^ t ^ . Initial film height 
uncertainties were, on average, ±5%. The dynamics of the film is recorded from t = ti 
onward. The angle of inclination was evaluated with image analysis and was evaluated 
within ±0.7 degrees. 

In a typical experiment, the downstream contact line advances in the direction of the 
flow (Fig.[ljd)). In the meantime, the film drains downstream and, in some cases, droplets 
were found to form and drip (Fig. [ljb)). However, for greater inclinations no droplets were 
observed to form (Fig. Wc)). Describing the full dynamics requires complex non-linear 
and time evolving tools^ and is beyond the scope of this article. Instead, we propose to 
simply count the number of droplets formed during an experiment, keeping track of their 
dripping time and location. In particular, we focus on the central part of the film and denote 
N the number of droplets dripping from this area. The region of interest, of area A^ is 
defined by scaling the perimeter inwards of a length 2n\/2£ c to avoid further complications 
at the contact line (see Supplementary Information for more details). To summarize, any 
experiment is fully characterized by the triplet: 

(hi/e c ,a,N) ex p (19) 

where N = N/Ai denotes the density of droplets in the center part of the flow. Those data 
points are reported in Fig. |IV| and are discussed in the next section. 
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As a side remark, note that the ratio hi/i c directly relates to the contact angle of the 
fluid on the substrate. In fact, the thickness of a partially wetting puddle sitting on an 
horizontal surface may be expressed as h p = 2^ c sin# c /2, in the limit of very long times 17 . 
By assuming the initial film thickness being the one of the puddle, the results may be 
represented as a function of 0 C instead of hi/i c . 


V. DISCUSSION 


As inferred from Fig. |6|a), the experiments divide into two categories: experiments that 
do not lead to any drop N = 0, and experiments during which several drops are formed 
N > 0. As anticipated, the largest drop-density experiments are the one with the smallest 
inclination of the substrate (a ~ 7r/2). In those cases, the flow is weak and does not affect 
the RTI. Increasing the inclinations (decreasing a when keeping other parameters fixed) 
yields a decrease in the number of drops. Interestingly, experiments reveal the existence of 
a critical angle a® xp under which no droplets are found. As seen from Fig. [6| this angle is 
a decreasing function of hi/i c . Those observations are now rationalized using the previous 
theoretical results. The absolute/convective transition angle a*(hi/£ c ) from eq. © is 
superimposed to the experimental results in Fig.[6ja) and Fig.Jf^b) in the form of a curved 
transparent surface and solid black line respectively. This theoretical prediction acts as a 
demarcation line between the two experimental regimes (N > 0 and N = 0) over the range 
of experimental parameters explored. 

On the one hand, this agreement was expected as the physical mechanism behind a 
convective instability is precisely that the flow outruns the instability, which cannot grow 
locally, and therefore drop formation is prevented. The geometry of the substrate, that is its 
inclination, stabilizes the flow. This type of geometrically induced stabilization is consistent 
with recent work on the dynamics of thin films on the underside of a cylinder where the 
curvature of the substrate, that is the local change in its slope, is shown to suppress the 
RTI if the film is thin enougtP^. 

On the other hand, we note the favorable agreement of theory with the experiments de¬ 
spite the significant differences between the experimentally obtained flow and the model 
depicted in Fig. [2j In particular, the model is one dimensional and our analysis is restricted 
to describing the linear behavior of the system. In reality, dripping is inherently non-linear 
and the experi ment is essentially three dimensional. Yet, as it often turns out in fluid 
mechanic d 1 2 3 4 5 6 7 8 9 10 11 12 * 19 *, the linear predictions prove to be an acceptable representation of the prob¬ 
lem far beyond their validity range, thereby helping in the rationalization of the observed 
physical reality. 
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